Superfluidity of two-dimensional excitons in flat and harmonic traps 
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Superfluid exciton density and superfluid transition (crossover) temperature are calculated for 2D 
excitons in large-size flat and harmonic traps. A generalized local density approximation for the 
Kosterlitz-Thouless theory is developed. 
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I. INTRODUCTION 

Two-dimensional (2D) indirect exciton systems in cou- 
pled quantum wells (CQWs) [l| (or in single QW in 
strong normal electric field [4I), in which electrons and 
holes are spatially separated, at low temperatures can 
possess condensation and superfluidity [J S JB| (or 
crystallization in some region of phase diagram [aS @])' 
jersistent electric currents and quasi- Josephson effects 
9|, as well as interesting coherent optical phenomena 
10| . There is a tremendous progress now in the observa- 
tion of collective exciton properties in CQWs (see reviews 

Due to the divergence of phase fluctuations, Bose con- 
densation of excitons (possible in 3D systems in ex- 
tended 2D exciton systems at T ^ is impossible p^ . 
But at low temperatures the system can possess quasi- 
long (power-law) off-diagonal order [l^ and superfluidity, 
which are destroyed by vortex pair dissociation and free 
vortex formation at Kosterlitz-Thouless (KT) tempera- 
ture [H, [1^ [131 • But Bose condensation is possible in 
finite 2D systems [2H. Besides, the trapping can increase 
the exciton density at fixed laser pumping, which leads 
to the growth of Bose condensation temperature. Several 
types of traps were proposed and experimentally realized: 

(i) " natural" traps created by (one-monolayer) well width 
fluctuations [l^, [H, [l^]; and regular, controlled traps 
such as (ii) laser-induced traps formed inside the annu- 
lus of laser pump [25] , (iii) traps formed by an inhomo- 
geneous deformation created by a needle stressing on the 
sample [26 , 27], (iv) mesa created by engineering width 
of layers or the barrier in CQWs, and (v) different kinds 
of electrostatic traps [H, [H, [29l . |30| . The interesting op- 
tical properties of exciton collective state are observed in 
luminescence measurements [H, H, [H [H HO, [H [H H 
(see also [ll!,[ll,[ll,[ll). 

In connection with the discussion above essential prob- 
lems arose concerning the nature of exciton collective 
state in 2D traps: (i) How the phase fluctuations and vor- 
tices manifest itself in sufficiently large but finite traps? 

(ii) How the type of the confining potential infiuences on 
the superfiuid density spatial distribution and its evolu- 
tion with the temperature? These problems the present 
paper are mainly addressed to. 

We describe quantitatively the superfluidity of 2D 



CQW excitons in 2D large-size flat trap by means of KT 
theory. For the investigation of the exciton superfluidity 
in 2D large-size harmonic trap we employ the local den- 
sity approximation [sS] and develop more correct gener- 
alized local density approximation for KT theory. 



II. SUPERFLUIDITY OF 2D EXCITONS IN A 
LARGE-SIZE HARMONIC TRAP: THE LOCAL 
DENSITY APPROXIMATION 

As the first step for study of 2D CQW exciton super- 
fluidity in a large-size harmonic trap we use the local den- 
sity approximation (LDA). LDA is valid if conflning po- 
tential of the trap varies sufficiently slowly on all micro- 
scopical characteristic scales of the problem, particularly, 
related to phonon and vortex excitations in the exciton 
system. Therefore, in LDA exciton superfluid properties 
near the trap point R coincides with properties of the 
homogeneous infinite superfluid with total density 
being equal to the (total) exciton density 71.4 (R) in the 
trap in the point R. 

A homogeneous infinite exciton system at fixed tem- 
perature T is a superfluid at high total densities but 
it is a normal fluid at low n^. As total exciton den- 
sity nj", when decreasing, goes through a critical one 
corresponding to exciton temperature T the super- 
fluid KT transition is occured in the system and super- 
fluid exciton density nf vanishes by the universal jump 
hf — 2mT/TTh^ [3^ (m is the exciton mass). 

So, the trapped excitons in LDA are superfluid (nor- 
mal) ones at the trap point R at which total exciton 
density nt(R) is greater (less) than the critical one n^, 
corresponding to exciton temperature T. Thus, in LDA 
at the sufficiently low temperatures at which the criti- 



cal total exciton density n"^ is smaller than maximal one 
n™ = nt{0) the exciton system in a symmetrical har- 
monic trap is divided into a superfluid circle R < Rs oi 
a radius Rs and normal ring R^ < R < L surrounding 
the circle (L is the Thomas-Fermi (TF) radius). In the 
superfluid circle the total exciton density nt{R) exceeds 
a critical value corresponding to exciton temperature 
T. But the superfluid density ns{R) is proved (see details 
in [19, 35]) to be higher than the universal jump n"^ for 
this temperature. In opposite, in the normal ring total 
exciton density nt{R) is smaller than and superfluid 
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density is equal to zero ns{R) = 0. In the circumference 
R = Rg which is the boundary line of the superfluid cir- 
cle and normal ring the local KT transition takes place 
on which the onset of pair vortex dissotiation and free 
vortex formation occurs. At the circumference R = R^ 
the superfluid density jumps from to zero. 

As the temperature grows critical total exciton density 
increases. So, the superfluid circle shrinkes and the 
normal ring thickens. At the temperature T^, at which 
the critical total density becomes equal to maximal total 
density in the whole trap n™ = nt(0), the superfluid cir- 
cle disappears and the trapped exciton system undergoes 
a qlobal transition to the normal state. 



III. SUPERFLUIDITY OF 2D EXCITONS IN A 
LARGE-SIZE FLAT TRAP 

In 2D homogeneous (flat) trap of sufficiently large size 
L the true, vortex-renormalized (VR), superfluid exci- 
ton density rig takes into account the renormalization 
("screening") by all vortex pairs in the system, i.e., by 
vortex pairs with the separation less than L. So, VR su- 
perfluid density Ug is expressed through "local", vortex- 
unrenormalized (VU), superfluid density ni as follows 



Us = ni/e{L), 

where e(r) is vortex pair dielectric fimction [l^ . 
In superfluid phase where there is a rare vortex 
the quantity e(L) according to KT model [l^ [3( 
the equation 



de{x) 
dx 



C(a)e''^"°-^o" 



(1) 

MM- 

Dair gas 
I obeys 



(2) 



with the boundary condition 

e(0) = 1. 



(3) 



Here x ~ ln(L/ro), rg is vortex core diameter (in CQWs 
To is of order of the distance between excitons [1^), 



2TTh^ni/mT, 



(4) 



and C(a) is a (nonuniversal) constant depending on mi- 
croscopical exciton superfluid properties; in 2D X-Y 
model it is equal to C(a) = Tr^ae"'^"/^. 

According to our numerical analysis of Eqs. H])-©, 
at low temperatures T (i.e., at large a; see Eq. (HI) VR 
superfluid density Ug is close to VU one n;. As T grows 
the ratio Us/ni decreases. At a critical temperature 
in sufficiently large traps (ln(L/ro) » 1) there is the 
crossover where rather sharp disappearance of the VR 
superfluid density takes place. The disappearance corre- 
sponds to the onset of pair vortex dissociation and free 
vortex formation [l^ , superfluid state being destroyed. 

Numerical integrating Eqs. ^ yield the following 
expression for the superfluid crossover temperature for 
CQW excitons in 2D flat trap of size L: 



T, 
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FIG. 1: The dependence of the superfluid crossover temper- 
ature of 2D GQW excitons in the flat trap on its size L 
{T^ = Tc\j^_^)- Thin solid line is the numerical intergation 
of the equation for vortex pair dielectric function e(L) in KT 
model, dash line is the fit e°°/e^ (according to Eq. ((6|), and 
dot line is the Kosterlitz's renormalization group calculation. 



where ni(T^) is the VU superfluid density of 2D excitons 
at the crossover (at T = T^) and the quantity is 



calculated numerically. The flt 



= ( 1 - 



(ln(L/ro) + A)^ 



(6) 



coincides excellently with the numerical calculation at 
large L (see below and Fig. 1); in Eq. ^ e°° is the criti- 
cal vortex pair dielectric constant in 2D inflnite superfluid 
at KT transition, h and A are nonuniversal quantities 
which, as i°° , depend on microscopical exciton system 
properties. In 2D X-Y model we obtain e°° = 1.134(9), 
A = 2.(93)-hA and b = 0.8(00). The quantity A = A(L) 
determines the correctness of KT model, so, within KT 
model one must assume that A = 0. The numerical cal- 
culation of the width of superfluid crossover region, in 
which vortex pairs should begin to dissociate, yields the 
following estimation ~ 1. 

From Eqs. (O, © we conclude that the approach used 
takes into account the second order on the logarithmically 
small quantity l/(ln(L/ro) + A). Thus, our calculation is 
correct in the logarithmic approximation (ln(L/ro) 3> 1), 
i.e., for sufficiently large traps. 

For the halfwidth of superfluid crossover region we ob- 
tain (see Eqs. (0 and (O) 



ATf 



(ln(L/ro) + A)3 



(7) 



Thus, the model adopted does not take into account ef- 
fects in the crossover region (because IST^ depends on 
the third order of ln(L/ro) ^ 1). In the crossover tem- 
perature region AT^^ the number of free vortices inside 
the trap area L"^ increases by e times, whereas at T = 
there is (in average) of order of one free vortex. 

Note that the result ([7]) for 2D differs essentially from 
the case of 3D Bose condensed systems where IST^ de- 
pends on a power of L. 
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FIG. 2: The dependence of the superfluid density of 2D CQW 
excitons in the flat trap of different size L on the ratio of 
exciton temperature T to the superfluid transition tempera- 
ture of 2D infinite systems. The regions of the superfluid 
crossovers on which the superfluid density universal jump is 
smeared out are shown by arrows. 



For the completeness of the calculation of VR super- 
fluid density Us and superfluid crossover temperature 
in the flat trap one must calculate the dependence 
of VU superfluid density ni on temperature T and to- 
tal density tit (see Eqs. (|T|), © and This cal- 
culation is realized in Refs. [38|, |39|. In this Letter 
we make a simplification. Namely, for spatially indi- 
rect CQW excitons the strong dipole-dipole repelling re- 
sults in a small thermal depletion of their VU super- 
fluid density (not the condensate) throughout the super- 
fluid phase (at T < T^) 38]. In a sufficiently large 
flat trap in which the condition > e°°/1.2 is hold 
(40| VU superfluid density ni thoughout the superfluid 
phase obeys to the inequality rit > ni > n* , where 
n* = ni{T^) = (0.93 0.98)nt for GaAs CQWs studied 
experimentally in Refs. [H [p^H, [H, M, 113, IH) 
and n* « 0.9nt in Refs. [Ijillil- So, ni is close to the 
total exciton density rit (for pure CQWs) [11]. Thus, for 
spatially indirect CQW excitons we can assume VU su- 
perfluid exciton density ni to be approximately equal to 
total one rit throughout the superfluid phase 



ni 



(8) 



This permits not only to simphfy Eqs. H]), (g]) and (O 



e{Ly 



2me 



mT 



T' 



(9) 



(10) 



but also to perform a complete calculation of the VR 
superfluid exciton density. (In Eqs. (O and pil)l the 
quantity Tq = 2'Kfi?nt/rn is the degeneration tempera- 
ture of the spin-polarized excitons; if the excitons are not 
spin-polarized and have gex spin degrees, the real exciton 
degeneration temperature is 2TrlT?nt/gex'm-, this quantity 



is in gex times less than Tq; gex = 4 for GaAs.) At typical 
total exciton density {in all spin degrees [4l|) nt = 2-10^*' 
cm^^ and exciton mass m — 0.22r7io (too is free electron 
mass) [2^ the indirect exciton superfluid crossover tem- 
perature in 2D flat trap is estimated as Te — 1.2 K for 
the trap size L — 5.6 /im, — 1.15 K for L = 43 /im, 
= 1.1 K for L = 14 cm and = 1.077 K for an infi- 
nite (L — 00) system (see Eqs. ^ and pUj) . degeneration 
temperature being equal to Tdeg = 1.26 K. 

In Fig. 2 we show the temperature dependence of the 
(VR) superfluid exciton density in 2D flat trap. The 
temperature unit is KT temperature of the 2D infinite 
superfiuid [11] 



'J-^L I 



L—oo 



(11) 



(here, as in Eq. (fTO]) . we assume ni{T^) — nt, see Eq. 
©). We see that the superfluid density universal jump 
fif = 2mT^ / nh^ = ntje^ is smeared out in a flat trap 
as compared to inflnite systems. Besides smearing, the 
confinement affects in growth of the superfluid crossover 
temperature. 

The dependence of superfluid crossover temperature 
Te on trap size L are depicted in Fig. 1. We see the fit ([6]) 
is much more better coincides with the exact KT-model 
numerical intergation of Eq. than the Kosterlitz's 
renormalization group calculation [20j 
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\n\L/ro) 



(12) 



(although, in the strict sence, Eqs. (O and ([T^ coincide 
in the logarithmic approximation). 



IV. SUPERFLUIDITY OF 2D EXCITONS IN A 
LARGE-SIZE HARMONIC TRAP: THE 
GENERALIZED LOCAL DENSITY 
APPROXIMATION 

The local density approximation for trapped excitons 
described in Sect. |TT]is valid only if the total exciton den- 
sity varies sufficiently slowly on both typical phonon and 
vortex scales of the trapped superfiuid. But the typical 
vortex scales of 2D homogeneous infinite superfiuid have 
the divergence near KT transition [20l. [42|. So, close to 
KT transition the vortex scales cannot be much smaller 
than the (finite) trap size (TF radius) L. Moreover, at 
sufficiently low temperatures the density of (thermally 
created) vortex pairs is negligible throughout the trap 
except the narrow region near the boundary of the sys- 
tem. Thus, for more accurate description of the exciton 
superfiuidity in 2D inhomogeneous traps the usual LDA 
must be modified. 

The generalization of KT rare vortex pair gas the- 
ory for inhomogeneous systems is not a simple prob- 
lem. Indeed, vortex charge q being proportional to square 
root y/ni{IVj of VU superfiuid exciton density rii(R) [l^ 
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varies as trap point R displaces. So, vortex pairs are no 
longer neutral dipoles. This result in an inapplicability, 
in the strict sense, of KT theory for a inhomogeneous 
case. However, this theory describes correctly the su- 
perfluid on the small length scales (much smaller than 
trap size L), on which the vortex charge varies slowly. 
On the scales of order of L the integration of Eq. ^ 
must be effectively cutted off. But as in a flat trap, 
in a harmonic trap vortex properties depend logarith- 
mically on the length scale. So, at sufficiently large L 
(ln(Z/ro) 3> 1) in the logarithmic approximation the ex- 
act value of the cutoff in integrating Eq. ^ is not im- 
potant. This results in appearing a possibility to inves- 
tigate quantitatively CQW exciton superfluid in a large- 
size harmonic trap with a logarithmically small error in 
the approximation which we will call generalized local 
density approximation (GLDA) for KT theory of weakly 
inhomogeneous systems. In GLDA, VR exciton super- 
fluid properties in 2D large-size harmonic trap near to 
trap point R are replaced by ones in 2D flat trap with 
total density equal to total density rit (R) in harmonic 
trap, flat trap size L being of order of harmonic trap size 

L. 

We obtain readily the exciton fluid in 2D symmetrical 
harmonic trap at low temperatures to divide in GLDA, 
as in LDA, into the superfluid circle 



nt(i?) > , UsiR) - ntiR) (R < R,) 
and normal ring 

nt{R)<fL^, n,(i?)wO {R^<R<L). 



(13) 



(14) 



Here R^ is the superfluid circle radius in the harmonic 
trap of size L, and (see Eq. (fTT 



2me 



-T 



(15) 



is (critical) total exciton density in 2D flat trap of 
size L ~ L, at which excitons undergo the superfluid 
crossover at temperature T. 

But near the circumference R Rs, where rit(R) w 
, which separates the superfluid circle and normal ring, 
there is a narrow region of the local KT crossover, on 
which the onset of the pair vortex dissociation and free 
vortex formation is gradual. This is one of distinctions 
between LDA and GLDA. 

The critical total density (fTS)) grows with the exciton 
temperature. So, superfluid circle radius R^ drops, the 
superfluid circle shrinks, and the normal ring becomes 
wider (see Eqs. (IT5|) - (IT5]) '). At the temperature (see Eq. 

m) 



VR superfluid component vanishes (in contrast to LDA) 
gradually. Above there is no superfluid circle, so the 
trapped excitons are in globally normal phase. 

We describe quantitatively the superfluidity of 2D 
CQW excitons in a large-size harmonic trap in GLDA 
for TF density profile [lllil 

nt{R) = ni{R) = n™(l - R^/L^)e{L - R), 

where 9{x) = if a; < and 6{x) = 1 if x > 0. In 
this case, we obtain in GLDA the following equations for 
superfluid circle radius R^, global superfluid crossover 
temperature , TF radius L, exciton chemical potential 
/I, maximal total exciton density n™ (in trap center), and 
exciton numbers in the superfluid circle Nsc and normal 
ring Nnr: 



hlLlo\fN't 



Nsc = Nt{l - {TlTtY)e{T^ - T), 



Nnr = Nt{{T/T^Ye{T^ -T) + 0{T - T^)), 



where ujo is the oscillator (trap) frequency, Xo — \Jhjmuio 
is the oscillator length, Nt is the total exciton num- 
ber in the trap (in all spin degrees [4l|), and C' = 
/im/27r?i^n™ = const is dimensionless exciton chemical 
potential, which we assume to be independent of T and 
n™. Local and global superfluid crossover halfwidths de- 
termining the correctness of the logarithmic approxima- 
tion (and hence, GLDA), as in a flat trap, have the third 
order on the logaritmically small quantity, i. e.. 



ARi 



2RLTL (ln(L/ro) + A)^ 



2me^ 



(16) 



AT' 



(ln(L/ro) + Ay 



at which the critical total exciton density is equal 
to the maximal one in whole the trap = nt(0) the 
excitons undergoes the global superfluid crossover where 



In Fig. 3 we show (VR) superfluid profiles for 2D har- 
monically trapped CQW excitons in GLDA and in LDA 
(see Fig. 3(a)) at different temperatures below (the 
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FIG. 3: (a) The vortex-renormalized superfluid profiles of 
2D CQW excitons in the harmonic trap of large radius L 
(ln(Z/ro) = 8) at different temperatures T in GLDA (thick 
lines) and in LDA (thin lines), (b) The vortex-unrenormalized 
superfl uid p rofiles in CQW structure (studied experimentally 



in Ref. []J,l221) for nT = 2.5 ■ 10^° cm"^ at the same temper- 
atures for corresponding lines as in (a), (c) Solid lines shows 
the temperature dependence of superfluid (thick) and super- 
fluid circle (thin) fractions. Dash lines denotes normal (thick) 
and normal ring (thin) fractions. 



temperature unit is the global superfluid transition tem- 



perature in LDA = f^\^^^ 
see Eqs. dH), dTO]), ([iT]) and VU superfluid pro- 

files a re p resented in Fig. 3(b) for CQW structure (see, 
6- g-, [l3[22|]) with exciton mass m = 0.22mo, dielectric 
constant e = 12.5, and effective electron- hole layer dis- 
tance d = 13.6 nm for trap center total exciton density 
being equal to = 2.5 • 10^° cm^^ One see the su- 
perfluid density on the superfluid circle boundary region 
to drop sharply far from the global superfluid transition 



(T 56 T|f ). Besides, there is little normal component in 
the superfluid circle at T 9^ T^^ (for a pure CQW struc- 
ture). The latter fact is visually depicted in Fig. 3(c), 
where the curve of the superfluid (normal) exciton num- 
ber in GLDA quite merge with that of exciton number 
in the superfluid circle (normal ring). 

Because the global superfluid crossover temperature 
Tjf for the harmonic trap of radius L coincides with the 
temperature for the flat trap of size L ^ L (see 
above), the dependence of the former on L coincides with 
that presented in Fig. 1 with the change T^^ 
yoo _^ foo^ .^^^ L-^L. 

V. CONCLUSION 



We have calculated the temperature dependence of su- 
perfluid density and the superfluid crossover temperature 
for 2D spatially indirect CQW excitons in a large-size flat 
trap as well as superfluid profiles and global superfluid 
crossover temperature for the excitons in a large-size har- 
monic trap. We reveal that harmonically trapped exci- 
ton system at low temperatures is divided into the su- 
perfluid circle (with superfluid component) and normal 
ring (without it). At their boundary there is the narrow 
crossover region at which the onset of vortex pair dissoci- 
ation and free vortex formation takes place. The super- 
fluid density goes rather sharply to zero in this region. As 
exciton temperature grows superfluid circle shrinkes to 
zero. This temperature corresponds to the global super- 
fluid transition (crossover) in the harmonic trap. At this 
temperature the Bose condensate also disappears (details 
will be published elsewhere). 
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